Fall and rise of small droplets on rough hydrophobic substrates 
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Liquid droplets on patterned hydrophobic substrates are typically observed either in the Wenzel 
or the Cassie state. Here we show that for droplets of comparable size to the roughness scale an 
additional local equilibrium state exists, where the droplet is immersed in the texture, but not yet 
contacts the bottom grooves. Upon evaporation, a droplet in this state enters the Cassie state, 
leading to a qualitatively new self-cleaning mechanism. The effect is of generic character and is 
expected to occur in any hydrophobic capillary wetting situation where a spherical liquid reservoir 
is involved. 

PACS numbers: 68. 08. Be, 68.35.Ct, 47.55. D- 
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The fact that roughness at the micrometer level can 
drastically increase the water-repellant properties of a 
hydrophobic substance has been known for a long time 
0, [H and is a frequent phenomenon in nature - the most 
prominent example being the Lotus leaf. Technical ad- 
vancement in surface fabrication and novel industrial ap- 
plications, such as self-cleaning materials, has brought 
the phenomenon of superhydrophobicity again in the fo- 
cus of research during the last years (see [3|, |4| for recent 
reviews). 

Until now, mostly droplets that are much larger than 
the typical roughness scale of the surface have been in- 
vestigated. In that case, the droplet either appears in the 
Wenzel state, where it completely wets the substrate 
or in the Cassie state, i.e. on top of the roughness struc- 
tures [2]. The fact that usually both the Cassie and the 
Wenzel state can be observed on the same substrate im- 
plies that both are separated by a free energy barrier, 
which can be overcome by external forces or kinetic en- 
ergy SH0. 

It is known, that the transition from the Cassie to the 
Wenzel state proceeds through the nucleation of contact 
between the liquid and the substrate at the grooves 
initiated e.g. via the increase of internal droplet pressure 
during evaporation 0, [Io|, EH- ^ * s unclear, however, 
how deep a droplet can really sink into the texture and 
how this process is modified for droplets of similar size as 
the roughness scale. Droplets of this length scale are not 
only important for a better understanding of the wetting 
properties of microscale systems [12| , but also play crucial 
roles in condensation processes on rough surfaces H, Q . 

The validity and applicability of the classical two-state 
picture of Wenzel and Cassie has been questioned several 
times [TBI, [l6|, [I?}- In particular, for the case of droplets 
that are of comparable size to the surface roughness - 
and which therefore directly feel the influence of the sur- 
face geometry - it can not be expected to hold a priori. 
Recent experimental studies, focusing explicitly on small 
droplets, confirmed at least the existence of the two equi- 
librium states 18, 
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FIG. 1: A small droplet on a hydrophobic pillar array, (a) 
shows the impaled state as a typical situation in the present 
simulations. In the analytical model (b), the droplet is as- 
sumed to be pinned at the edges of the pillars, i.e. it has a 
fixed base radius a. For a given droplet size, the only way to 
minimize the overall free energy is thus through a change in 
the penetration depth p. This also fixes the apparent contact 
angle 0. 



Here we show, both analytically and via lattice Boltz- 
mann (LB) simulations, that for a droplet of comparable 
size to the roughness scale, there exists a further generic 
state (hereafter called impaled state), characterized by 
almost complete immersion of the droplet into the tex- 
ture, yet not touching the bottom of the grooves. To 
the best of our knowledge, such a state has not been re- 
ported before. It is important to realize that this state is 
different from the "partially impaled" conformations of 
the Cassie state, where the liquid- vapor interface below 
the macroscopic contact line is curved, but the droplet 
essentially remains on top of the asperities [Io|, [III • This 
statement is supported by the fact that we indeed find 
a coexistence regime for the new impaled state and the 
(partially impaled) Cassie state. Remarkably, by virtue 
of this new state, a droplet can reenter a Cassie state 
upon, e.g., evaporation. 

In this work, we choose the simplest geometry that can 
model the roughness of a surface, namely, a regular array 
of cuboidal pillars with width 6, height h and spacing d 
(Fig.[TJ). The intrinsic hydrophobicity of the fiat parts of 
the surface is described by the Young contact angle Oy. 
In this work we only focus on the case Oy > 90°. Gravity 
will be neglected throughout, as we only consider droplets 
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that are smaller than the capillary length (<~2.7 mm for 
water). 

In contrast to standard single phase LB models [2]], 
I22I [23[, we employ here a free energy based two-phase 
(liquid- vapor) LB approach 24, 25[. Details of the algo- 
rithm and the parameters used can be found i n 12 61 127[ . 
The relation between LB and physical units [26[ shows 
that our simulated droplets would be of micron scale. 
However, since it can be argued that the actual value 
of the physical viscosity is not important in the present 
case, it can be used to tune the unit of length, allowing 
the simulation results to be applied to a broad range of 
length scales, as long as thermal fluctuations and gravity 
can be disregarded (see also below). 

In the analytical model (Fig. [TJd), we assume the part 
of the droplet located above the pillars to be a spheri- 
cal cap with base a = b/2 + d = RsinO (with R being 
the radius of the cap and the apparent contact angle). 
The impaled part is approximated as a hollow cylindrical 
column with radius a and height p (penetration depth). 
The macroscopic contact line of the droplet is assumed 
to remain pinned at the outer pillars. Note that in this 
model, the Wenzel state would correspond to p = h and 
the Cassie state to p = 0. Therefore, it also neglects a 
possible finite penetration depth of a droplet in a "par- 
tially impaled" Cassie state. As further simulations Q 
have shown, the symmetric droplet configuration consid- 
ered here is stable against moderate perturbations. 

The total droplet volume is fixed, Vtot = 
const = V sph (0) + Vcy\(p), with V sph = |7ra 3 (2 - 
3 cos + cos 3 0) I sin 3 the volume of the cap and 
V cy \ = (na 2 — b 2 )p the volume of the penetrating cylin- 
der. In the following, we will refer to the effective droplet 
radius i? e ff (47r/3i? 3 ff = Vtot) as the appropriate measure 
of droplet size, p is the free variable and 6 is determined 
via the fixed volume condition. 

Since the volume of the drop (and the temperature) 
is constant, only surface energy contributions play a role 
for a change in the total free energy. The free energy 
f(p) of the model droplet, neglecting gravity and terms 
associated with the Wenzel transition, and normalizing 
such that /(0) = 0, is then given by 

f(p) = ^lv S sph (p) - S sph (0) + 5 cy i,Lv(p) - 8bp cos Y 

(1) 

Here, 5 sp h = 27ra 2 (l — cos#)/ sin 2 is the surface area 
of the spherical cap, £ C yi,LV = (2™ — 4b)p is the 
lateral liquid-vapor surface area of the cylinder and 
— au/Sbpcosdy is the (positive, since 0y > 90°) free en- 
ergy associated with the wetting of the (eight) side walls 
of the pillars. 

Figures [2^i,b show the dependence of the free-energy 
on the penetration depth p for varying Young contact 
angles and droplet sizes. Several interesting observations 
can be made: Firstly, as also found in the case of droplets 
large compared to the roughness scale d, [Io|, EH, , the 



stability of the Cassie state, determined by the slope of 
/ at p = 0, depends not only on the contact angle but 
also on the size of the droplet. 

The novel feature is the appearance of a local mini- 
mum of the free energy at large penetration depths, ex- 
isting in addition to the possible minimum associated 
with the Cassie (and Wenzel) state. From the condition 
df I dp = 0, which is easily evaluated with the help of the 
fixed volume constraint, a necessary condition for the ex- 
istence of a minimum of the free energy arises, namely, 
0y < arccos( — \ + ^). Note that this relation does not 
depend on the droplet size. 

The origin of this new state can be easily understood 
by imagining the pillars to represent a (partly open) hy- 
drophobic capillary tube, wetted by a small droplet that 
is placed at its entry. In this situation, the local equilib- 
rium state of the droplet is a consequence of the balance 
between the Laplace pressure within the spherical cap 
(pushing the droplet into the capillary) and an opposing 
capillary force due to the hydrophobic substrate. 

Indeed, splitting up the free energy (Eq. [T]) into the 
contributions of the spherical cap and the remaining 
"capillary" part (Fig. [St), shows that the local minimum 
of / is ultimately a consequence of the non-linear de- 
pendence of the spherical cap surface on p. This simple 
argument suggests, that the intermediate minimum con- 
stitutes a generic state of a droplet, occurring in any 
situation of filling hydrophobic capillaries by a spherical 
liquid reservoir. Indeed, further simulations using vari- 
ous surface geometries (e.g. omitting the central pillar) 
clearly underline this assertion [2^ . 

Fig. [3] provides a test of our analytical free energy 
curves, thereby confirming also the existence of the in- 
termediate minimum. As predicted by the corresponding 
theoretical free energy curve (Fig. EK), the equilibrium 
position of a droplet depends on where it is placed at the 
beginning: A droplet deposited close to the top of the 
pillar array moves further to the top (Fig. 03d), while a 
droplet that initially penetrates deeper into the grooves 
becomes trapped in the intermediate minimum of the free 
energy (Figs. [3fc,d). 

Figure [4] depicts a phase diagram displaying the theo- 
retically expected regions of existence for the Cassie and 
the impaled state. The Cassie state is stable for the set 
of points (0y,i2eff) that fulfill df/dp\ p= o > 0, while the 
phase boundary for the impaled state is determined from 
df /dp = in combination with d 2 f /dp 2 > 0. 

Interestingly, the stability region of the Cassie state 
shows a characteristic shape, which is also largely inde- 
pendent of the geometry (Fig. Sj): There exists a certain 
droplet size where the Cassie state is unstable for a max- 
imal range in contact angle, and both towards larger and 
smaller radii the stability region increases. 

Noting that in the phase diagram a quasi-statically 
(2^ evaporating droplet would move on a vertical line 
from large towards small i? e ff, we infer the existence of 
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FIG. 2: Predictions of the analytical model, (a) Dependence of the free energy //47ra 2 aLv on the penetration depth p and 
the Young contact angle Oy for a droplet of fixed size i? e ff = a. (The inset shows a magnification of the curve for Oy = 108°.) 
(b) Dependence of the free energy //47ra 2 aLv on the penetration depth p and the droplet size i? e ff for a fixed contact angle 
of Oy — 100°. (c) Contributions to the free energy / / 'A.na 2 o"lv for i? e ff = a and Oy — 100°. / sp h is the surface free energy of 
the liquid- vapor interface of the spherical cap, f cap is the free energy due to the wetting of the "capillary" constituted by the 
pillars, and / = / sp h + /cap is the total free energy. All curves in (a-c) are given for b/d = 1 and plotted up to a value of p 
where no further volume is left in the spherical cap. 




FIG. 3: Droplet behavior for different initial configurations. We take a fixed droplet size of i? e ff = a, a contact angle of 
0y = 106° and b = d = 12 LB units, (a) shows the theoretical free energy //47ra 2 aLv in dependence of the penetration depth 
p. In (b-d), the results of the LB simulations are displayed. The labeled dots in (a) mark the initial penetration depth of the 
droplet in the cases (b-d). The arrows indicate the observed time evolution from the initial state to the final configuration. 
The predicted new impaled equilibrium state is indeed observed in the panels (c,d). 
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FIG. 4: Phase diagram: Regions of stability for the Cassie 
and the impaled state (intermediate minimum) as predicted 
by the analytical model. Focusing on the case of d/b = 1, 
the impaled state is expected to be stable in the complete 
shaded region, while the Cassie state is predicted to be stable 
right to the dashed curve. Dotted curves give the theoretical 
predictions for other values of d/b. The arrow illustrates the 
path of a quasi-statically evaporating droplet as it enters the 
region of an unstable Cassie state (o), and finally goes over 
to a stable Cassie state again (•). 



a reentrant transition: A droplet, initially existing in a 
stable Cassie state, can (following the arrow in Fig.HJ) be- 



come unstable upon a reduction of its size and thereupon 
adopt an impaled state (o). However, due to the fact that 
the position of the local free energy minimum shifts to- 
wards smaller p with decreasing volume (Fig.[2j3), further 
evaporation will always result in the droplet to re-appear 
in a stable Cassie state (•). 

This remarkable prediction is indeed confirmed by sim- 
ulation (Figure In the beginning (Fig.[5fci), the droplet 
resides in a ("partially impaled") Cassie state and only 
slightly increases its penetration depth (thus confirming 
[id |20[). However, after reaching a critical size, the 
droplet suddenly penetrates into the pillar grooves and 
goes over to the intermediate minimum of the free en- 
ergy (Fig. [SJd). During the impaled phase (Fig. [5)3, c) 
the droplet gradually climbs up the pillars again, still re- 
siding in the local minimum. Note that its penetration 
depth is in nice agreement with the analytical predictions 

(Fig. BO- 

According to the common understanding of self- 
cleaning, impalement is considered unfavorable and the 
droplet cleans the surface by rolling over the top of the 
texture. Interestingly, the reentrant transition provides 
a qualitatively new self-cleaning mechanism, since it al- 
lows to clean not only the top of the substrate, but also 
its inner parts. 




FIG. 5: Reentrant transition through evaporation of a droplet. In the simulation, the droplet is initially (t=0) prepared in a 
stable ("partially impaled") Cassie state (a). In the course of the quasi-static evaporation process, the droplet becomes unstable 
and adopts the impaled state (b), from where it gradually climbs up the pillars (c) until it reaches again a Cassie state (d). (e) 
shows the penetration depth of the lower droplet interface as it is observed in the simulation (•) and predicted by the analytical 
model ( — ). Simulation parameters: contact angle By = 103°, initial droplet size R e ff = 1.2a, b = d = 12 LB units. Evaporation 
is simulated by reducing the mass of the vapor phase at a rate of 5 x 10 _7 /LB time steps across a xy— plane close to the top 
of the simulation box. 



The existence of a reentrant transition can also explain 
some recent experimental observations, that found that 
small evaporating droplets indeed tend to remain close 
to the top of the substrate [l, [ll[ and not get trapped 
inside of the texture. 

It is important to stress the scale invariance of Eq. [TJ 
where thermal fluctuations and gravity are neglected. 
This may, firstly, significantly widen the range of possi- 
ble applications, and secondly, considerably simplify ex- 
perimental verification of our predictions. One could e.g. 
study sub-millimetric drops on pillars of comparable size, 
thereby avoiding problems such as preparation and fast 
evaporation of micron sized droplets. 

In conclusion, via analytical calculations and numerical 
simulations, we have uncovered the existence of an addi- 
tional local equilibrium state, where the droplet partially 
wets the inside of the grooves, yet is not touching the 
bottom of the substrate. This finding qualitatively mod- 
ifies the two-state paradigm of Wenzel and Cassie states 
in the case of droplets of comparable size to the surface 
roughness. Interestingly, droplets in this new state ap- 
pear to have Wenzel-like properties (e.g. small apparent 
contact angles) - but on the other hand possess the in- 
herent capability to reenter the Cassie state again. We 
have demonstrated that our results are largely indepen- 
dent of the particular surface geometry and are expected 
to hold whenever a hydrophobic capillary is wetted by a 
small droplet. 
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